AO-A047  •!«  AIR  FORCE  ZNST  OF  TECH  WRIOHT-PATTERSON  AFB  OHIO  F/«  20/4 

aircraft  compensator  OESI0N  METHODS. (U) 

JUN  77  T U RIMS 


UNCLASSIFIED AFIT-CI-77-i5  NL 


AO 

A047ei4 


; iw  iiws 

flOC  FILE  COPY  ADA047814 


Tom  L.  Riggs,  Jr. 


Certificate  of  Approval: 


Electrical 


Professor 
Engineering 


S.  Boland,  I!!,  C^airrar. 


Associate  Professor 
Electrical  Engineering 


D.  W.  Russell,  Professor 
Electrical  Engineering 


"Paul  F.  Parks,  Dean 
Graduate  School 


Tom  L.  Riggs,  Jr. 


A Thesis 
Submitted  to 
the  Graduate  Faculty  of 
Auburn  University 
in  Partial  Fulfillment  of  the 
Requirements  for  the 
Degree  of 

Master  of  Electrical  Engineering 


VITA 


Tom  Lance  Riggs,  Jr.,  son  of  Tom  Lance  and  Fern  Marie  (Lowry) 
Riggs,  was  born  in  Washington,  D.  C.  on  April  21,  1947.  He  attended 
public  school  in  Washington,  D.  C.  until  the  age  of  15  when  he  and  his 
family  moved  to  Arlington,  Virginia.  He  completed  his  secondary  edu- 
cation in  Arlington  County  Public  Schools.  In  April,  1966,  he  enlisted 
in  the  U.  S.  Air  Force.  He  served  in  the  U.  S.  Air  Force  as  an  elec- 
tronics technician  until  1973  whereupon  he  was  selected  by  the  U.  S. 

Air  Force  to  attend  Auburn  University  to  attain  a baccalaureate  degree 
in  Electrical  Engineering  with  the  ultimate  goal  of  being  commissioned 
a 2nd  Lieutenant  in  the  U.  S.  Air  Force.  He  entered  Auburn  University 
in  September,  1973  and  received  the  degree  of  Bachelor  of  Electrical 
Engineering  in  March,  1976.  By  delaying  his  commission  for  one  year, 
the  U.  S.  Air  Force  selected  him  to  enter  the  Graduate  School,  Auburn 
University,  in  March,  1976.  He  married  Patricia  Frances,  daughter  of 
John  William  and  Margaret  May  (Donnison)  Davis  of  Manchester,  England 
in  December,  1969.  They  have  three  children.  Dawn  Michelle,  Karen 
Nichole,  and  Kimberly  Annette. 


iv 


r 


THESIS  ABSTRACT 

AIRCRAFT  COMPENSATOR  DESIGN  METHODS 
Tom  L.  Riggs,  Jr. 

Master  of  Electrical  Engineering,  June  7,  1977 
(B.E.E.,  Auburn  University,  1976) 

132  Typed  Pages 

Directed  by  Joseph  S.  Boland,  III 

> 

^ In  this  paper  methods  using  linear  analysis  are  presented  for  the 
design  of  longitudinal  flight  compensators.  The  compensators  employ 
state  feedback  to  force  the  aircraft  to  respond  in  the  desired  manner. 

The  design  process  involves  analysis  of  the  uncompensated  aircraft, 
construction  of  a model  which  has  the  desired  response,  and  two  algo- 
rithms for  designing  the  compensator. 

A method  is  presented  for  contriving  the  desired  model  from  handling 
qualities  performance  criteria.  The  method  is  straightforward  and 
results  in  the  exact  desired  short  period  response,'  however,  the  phugoid 
(long  period)  response  is  less  predictable. 

The  compensator  design  algorithms  are  easily  implemented  into 
computer  programs.  The  algorithms  require  a minimum  of  human/computer 
interaction  and  solutions  are  assured  for  controllable  systems. 

Compensators  for  two  aircraft  are  designed  to  show  the  techniques 
presented.  In  both  examples  the  compensated  system  response  is  exactly 
the  same  as  the  model  response.  ‘ 


L 


/ V 


TABLE  OF  CONTENTS 


LIST  OF  TABLES vii 

LIST  OF  FIGURES viii 

I.  INTRODUCTION 1 

II.  DEVELOPMENT  OF  AIRCRAFT  LONGITUDINAL  EQUATIONS  AND 

PERFORMANCE  CRITERIA 5 

Development  of  Equations 
Short  Period  Approximation 
Phugoid  Approximation 
Performance  Criteria 
Algorithm  for  Generating  Model 

III.  DEVELOPMENT  OF  COMPENSATOR  ALGORITHMS  31 

Concept  of  Equi control  lability 
Curran’s  Algorithm 
Alternate  Method 

IV.  APPLICATION  OF  COMPENSATOR  DESIGN  ALGORITHMS  TO  THE 

AIRCRAFT  LONGITUDINAL  CONTROL  PROBLEM  76 

V.  SUMMARY  AND  CONCLUSIONS 119 

BIBLIOGRAPHY 122 


LIST  OF  TABLES 


-1  Dimensional  and  Nondimensional  Aerodynamic  Forces 


vii 


LIST  OF  FIGURES 


2-1  Angles  of  Longitudinal  Flight 9 

2-2  The  A-Matrix  for  the  Aerodynamic  Equations  of  Longitudinal 

Flight  in  State  Space  Form 12 

2-3  The  B-Matrix  for  the  Aerodynamic  Equations  of  Longitudinal 

Flight  in  State  Space  Form 13 

2-4  Performance  Criteria  Plot  When  < 15 23 

2- 5  Performance  Criteria  Plot  When  > 15 23 

3- 1  System  Block  Diagram  Representation  of  Mathematic  Aircraft 

Model 32 

3-2  System  Block  Diagram  of  Proposed  Form  of  Compensation 32 

3-3  Block  Diagram  Representation  of  State  Equations  for  the 

Equicontrollable  System  Given  in  Example  3.1  37 

3-4  Block  Diagram  Representation  of  State  Equations  for  the 

System  in  Example  3.2 39 

3-5  Block  Diagram  Showing  the  Method  of  Connecting  the  Added 

States  to  the  Driginal  System 45 

3-6  Block  Diagram  Representation  of  the  Compensator 50 

3-7  Block  Diagram  Representation  of  the  State  Equations  for  the 

System  in  Canonical  Form  in  Example  3.3 54 

3-8  Final  Form  of  Compensated  System  in  Example  3.3 60 

3- 9  Block  Diagram  for  Compensated  System  Given  in  Proposition  3.  . 75 

4- 1  Handling  Qualities  for  System  in  Example  4.1  81 

4-2  Change  in  Velocity  vs.  Time  for  System  in  Example  4.1 90 

4-3  Change  in  Angle  of  Attack  vs.  Time  for  System  in 

Example  4.1 

viii 


I 

4-4  Change  n Pitch  Angle  vs.  Time  For  System  in  Example  4.1.  . . 92 

4-5  Pitch  Rate  vs.  Time  for  System  in  Example  4.1 93 

4-6  Change  in  Altitude  vs.  Time  For  System  in  Example  4.1  ...  . 94 

4-7  Change  in  Velocity  vs.  Time  For  System  in  Example  4.1  ...  . 95 

4-8  Change  in  Angle  of  Attack  vs.  Time  For  System  in 

Example  4.1 96 

4-9  Change  in  Pitch  Angle  vs.  Time  For  System  in  Example  4.1.  . . 97 

4-10  Pitch  Rate  vs.  Time  For  System  in  Example  4.1 98 

4-11  Change  in  Altitude  vs.  Time  For  System  in  Example  4.1  ....  99 

I 

4-12  Handling  Qualities  For  System  in  Example  4.2 102  i 

4-13  Change  in  Velocity  vs.  Time  For  System  in  Example  4.2  ....  105 

4-14  Change  in  Angle  of  Attack  vs.  Time  For  System  in 

Example  4.2 106 

4-15  Change  in  Pitch  Angle  vs.  Time  For  System  in  Example  4.2.  . . 107 

4-16  Pitch  Rate  vs.  Time  For  System  in  Example  4.2 108 

4-17  Change  in  Altitude  vs.  Time  For  System  in  Example  4.2  ....  109  ^ 

4-18  Change  in  Velocity  vs.  Time  For  System  in  Example  4.2  ....  110  | 

4-19  Change  in  Angle  of  Attack  vs.  Time  For  System  in 

Example  4.2 Ill 

4-20  Change  in  Pitch  Angle  vs.  Time  For  System  in  Example  4.2.  . . 112 

4-21  Pitch  Rate  vs.  Time  for  System  in  Example  4.2 113 

4-22  Change  in  Altitude  vs.  Time  for  System  in  Example  4.2  ...  . 114 

4-23  Handling  Qualities  of  Compensated  System  in  Example  4.2  For 
Velocity  = 330  fps  and  Altitude  Ranging  from  Sea  Level  to 
30,000  ft 115 

4-24  Handling  Qualities  of  Compensated  System  in  Example  4.2  For 
Velocity  = 880  fps  and  Altitude  Ranging  from  Sea  Level  to 
30,000  ft 115 


ix 


4-25  Handling  Qualities  of  Compensated  System  in  Example  4.2  For 
Velocity  = 440  fps  and  Altitude  Ranging  From  Sea  Level  to 
30,000  ft 116 


Velocity  = 440  fps  and  Altitude  Ranging  From  Sea  Level  to 

30,000  ft 116 

4-26  Handling  Qualities  of  Compensated  System  in  Example  4.2  For 
Velocity  = 770  fps  and  Altitude  Ranging  From  Sea  Level  to 

30,000  ft 116 


4-27  Handling  Qualities  for  Compensated  System  in  Example  4.2  For 
Velocity  = 550  fps  and  Altitude  Ranging  From  10,000  to 

30.000  ft 117 

4-28  Handling  Qualities  for  Compensated  System  in  Example  4.2  For 
Velocity  = 550  fps  and  Altitude  Ranging  From  Sea  Level  to 

5.000  ft 117 

4-29  Handling  Qualities  for  Compensated  System  in  Example  4.2  For 
Velocity  = 660  fps  and  Altitude  Ranging  From  20,000  to 

30.000  ft 118 

4-30  Handling  Qualities  for  Compensated  System  in  Example  4.2  For 
Velocity  = 660  fps  and  Altitude  Ranging  From  Sea  Level  to 

15.000  ft 118 


X 


1 

, 

'ti 

I.  INTRODUCTION 

The  modern  aircraft  is  capable  of  flight  over  a wide  range  of 
velocities  and  altitudes.  This  increased  capability  has  resulted  in  a 
deterioration  in  airplane  stability  manifested  by  an  increase  in  the 
airframe  natural  frequencies  and  a decrease  in  the  airframe  damping  [1]. 

In  order  to  improve  the  stability  of  the  modern  aircraft,  automatic 
control  devices  have  been  devised.  These  control  devices,  often  called 
Stability  Augmentation  Systems  (SAS),  operate  almost  universally  by 
sensing  the  airframe  motions  and  then  moving  a control  surface  to 

( 

oppose  the  airframe  motion. 

The  problem  of  designing  an  SAS  is  complex.  The  equations  which 
represent  the  mathematical  model  of  the  airframe  are  nonlinear  differ- 
ential equations.  For  reasons  of  comfort,  safety  and  mission  perform- 
ance, it  is  desirable  that  the  airplane  fly  along  a smooth  path.  It 
is  reasonable  to  consider  departures  from  the  path  as  small  perturba- 
tions. This  assumption  permits  considerable  simplification  of  the 
airframe  equations  of  motion  by  reducing  the  model  to  two  independent 
linear  differential  systems,  one  representing  longitudinal  flight  and 
the  other  representing  lateral  flight. 

The  design  problem  has  now  been  broken  down  into  two  independent 
problems,  one  for  longitudinal  flight  and  the  other  for  lateral  flight. 

Although  the  problem  has  been  simplified  greatly,  the  design  of  the 
SAS  for  the  two  flight  configurations  is  by  no  means  trivial.  For 

i; 
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instance,  the  small  disturbance  mathematical  model  for  longitudinal 
flight  is  a 5th  order  linear  differential  system  with  multiple  inputs 
and  outputs.  In  order  to  design  an  SAS  for  longitudinal  flight,  the 
system  open  loop  transfer  function  and  the  desired  transfer  function 
must  be  known.  Then  using  some  algorithm  the  compensator  is  designed 
such  that  the  compensated  system  transfer  function  is  the  same  as  the 
desired  transfer  function. 

One  method  for  designing  the  SAS  is  by  determining  the  transfer 
function  for  a particular  state  with  respect  to  some  input  and  design- 
ing a compensator  that  places  the  poles  and  zeroes  of  the  transfer 
function  at  some  desired  location.  However  this  method  is  impractical 
for  even  moderate  size  systems  because  of  coupling  between  states. 

A more  practical  method  is  to  generate  a model,  normally  of  the 
same  dimension  as  the  plant  and  often  of  the  same  structure  as  the 
plant,  that  has  the  overall  desired  response  for  the  system.  The 
problem  now  reduces  to  determining  a compensator  that  will  force  the 
states  of  the  plant  to  equal  the  states  of  the  model.  The  bulk  of  the 
work  in  attempting  to  solve  this  problem  has  used  optimal  control 
theory.  To  do  this  requires  the  generation  of  a cost  function.  Most 
papers  written  in  this  area  specify  a quadratic  performance  measure 
as  given  in  Equation  (1-1). 

j .y  [(Xp  - 

0 


(1-1) 
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where 

= states  of  the  plant 
P 

x„  = states  of  the  model 
m 

u = control  applied 

* 

and  Q and  R are  nxn  and  mxm  weighting  matrices  respectively.  This 
problem  has  been  solved  by  Ryanski,  Reynolds,  and  Shed  [2]  and  Windsor 
and  Roy  [3]  among  others.  However  the  key  difficulty  in  all  of  these 
solutions  is  the  lack  of  definition  of  the  matrices  Q and  R. 

R.  T.  Curran  [4]  developed  an  algorithmic  approach  to  the  design 
of  the  compensator.  By  performing  transformations  on  the  plant  and 
model  systems  he  transforms  the  system  to  a special  canonical  form 
which  allows  him  to  easily  determine  the  compensator  which  equates  the 
transfer  functions  of  the  model  and  compensated  plant.  Since  the 
transfer  function  for  a given  system  is  not  dependent  on  the  coordinate 
system  of  the  state  equations,  a solution  to  the  problem  is  assured. 
Once  the  compensator  is  determined  it  is  transformed  back  to  the 
original  coordinate  system. 

The  aircraft  control  problem  lends  itself  to  Curran’s  method 
because  of  the  form  of  the  aircraft  system  matrices.  In  this  paper, 
Curran's  method  is  fully  developed  in  relation  to  the  aircraft  longi- 
tudinal control  problem.  Curran's  algorithm  is  simplified  and  extended 
to  include  a wider  range  of  systems.  This  paper  shows  that  a solution 
to  the  aircraft  longitudinal  control  problem  via  Curran's  method  is 


II.  DEVELOPMENT  OF  AIRCRAFT  LONGITUDINAL  EQUATIONS 
AND  PERFORMANCE  CRITERIA 


Aerodynamic  forces  are,  in  general,  roughly  proportional  to 

p p 

pV  t where  p is  the  density  of  air,  V is  the  velocity  relative  to  the 
local  environment,  and  «,  is  a characteristic  linear  dimension.  It  is 
therefore  universal  practice  to  obtain  nondimensional  force  coefficients 
by  dividing  the  aerodynamic  forces  by  a factor  proportional  to  the 
above  quantity.  For  example,  the  weight  coefficient  is  given  by 
Cy^  = 2mg/pV^S  where  mg  is  the  weight  of  the  aircraft  and  S is  usually 
the  wing  surface  area.  Table  2-1  lists  the  aerodynamic  forces  related 
to  longitudinal  flight,  the  small  disturbance  divisor  and  the  resul- 
tant nondimensional  quantity  [6].  The  constants  in  Table  2-1  are 
defined  below. 

Vg  = constant  reference  velocity 
Pg  = constant  reference  density 

c = longitudinal  reference  length  (wing  mean  chord) 

S = reference  wing  plan  area 

Using  the  nondimensional  quantities  of  Table  2-1,  a set  of  1st 
order  linear  nondimensional  dynamic  equations  can  be  written  which  will 
describe  the  aircraft's  longitudinal  motion.  The  equations  were  taken 
from  Etkin's  Dynamics  of  Atmospheric  Flight  [6]  and  are  given  by 
Equations  (2-1)  through  (2-5). 
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Table  2-1.  Dimensional  and  Nondimensional 
Aerodynamic  Forces. 


Dimensional 

Quantity 

Symbol 

Small  Disturbance 

Divisor 

Nondimensional 

Quantity 

Thrust,  Drag,  Lift 

O 

1— 

JspV^S 

Weight 

mg 

hpyh 

Cw 

Moment 

M 

hp^hc 

Cm 

Pitch  Rate 

q 

2Vg/c 

q 

Mass 

m 

hp^Sc 

V 

Inertia 

’y 

PgS(F/2)^ 

■y 

Density 

P 

Pe 

P 

Velocity 

V 

''e 

V 

Altitude 

z 

c/2 

z 
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2ii(DV)  = [C^^cos(a^)  - Cp^  + 2C^  sin(Yg)]AV 

- [Cp  + C,  sin(a^)]Aa  - [C^  COSYg]AY 
a e • e 

+ [C,  cos(o-r)  - sin(Y  )]aZ 

'z  9z  ”e  ® 

+ AC-j-  cos(ay)  - ACp 


(2-1) 


•(2u+CLj(Da)  = [C-p  sintay)  + + 2C^  cos(Yg)]AV 

'"a  V V e 

+ [C.  + COs(a_)]Aa  + [C.  - 2ii]q 

'e 

+ [C^  sin(Yg)]AY 

+ [C,  + C,  sin(aT)  + ^ C cos(y_)]aZ 

'z  ' 8z  ”e  ® 

+ aC,  sin{a^)  + aC, 

'c  c 


(2-2) 


I Dq  - C Da  = C aV  + C Aa  + C q + C aZ  + AC  (2-3) 
y ^ m. 


Da  + Dy  = q 


(2-4) 


DZ  = - sin(Yg)AV  - cos(Yg)AY  - SiOYg 


(2-5) 


The  symbols  in  the  Equations  (2-1)  through  (2-5)  are  defined  below. 


e = reference  state 


De  = ^ e 

2Vg  dt  2Vg 


(2-6) 
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aC  = nondimensional.ized  command  input  e 

a-  = angle  between  the  thrust  vector  and  the  longitudinal 
body  axis 

Yg  = steady  state  angle  of  climb 
Often  to  simplify  the  equations,  a new  angle,  0,  is  defined  such 

that 


ct> 

11 

Q 

+ 

-< 

(2-7) 

Figure  2-1  illustrates  the  angles  of  longitudinal  flight.  For  simu- 
lation purposes  it  is  assumed  that  at  time  t=0  the  aircraft  is  in 

equilibrium  and  is  flying  straight  and  level.  Equilibrium  of  the 

reference  state  connects  the  thrust,  drag,  lift,  and  weight  coefficients 

by  the  following  relations 

Ct^cosot  - = C^^sinyg 

(2-8) 

Cf^sincxT  + Cl^  = C^^cosYe 

(2-9) 

Level  flight  is  defined  as 

II 

O 

(2-10) 

If  the  thrust,  drag,  lift,  and  moment  all  vary  with  the  air  density 
when  the  speed  is  constant  then,  it  is  a good  assumption  [6]  that 


(2-n) 


Substituting  Equations  (2-6)  through  (2-11)  and  the  definitions  for 
V,  q,  and  Z into  Equations  (2-1)  through  (2-5),  the  longitudinal 
equations  become 

. V V ^ 

AV  = — 3 (Ct  COSaj  - Crj  )aV  + — ^ (Cl  - Cn  )Aa 
pc  'v  ' pc  ‘-e 


y 2 y 2 

+ — ^ (ACt  COSa-r  “ACpj  ) 

pc  "e  pc  'c  ' ‘^c 


(2-12) 


, Ct  sina_  + Cl  +20.  2Vp  Ci  + Cn 

= - I (-"^  T . Lv 

c 2p  + Cl_  c 2p  + Cl_ 

a a 


- Cl  1 /4Veyaa\r  c^^ 

2p  + C,  'W  Aaz/  2p  + C,  . 

dJ  L *-0. 


^'^e  ^Cj^sina-j-  + aC^^ 
c 2p  + C^, 


(2-13) 


0 = q 


(2-14) 


pt  ^ 

Lc J i 


(Ct  sinaj  + Cl  + 2Cv^  ) 

.C  - C — aV 

y '"v  ’’’d  2p  + Cl. 


2Ve  ' 1 r (^L  + Cp  ) ] 

. ^ '"a  ‘‘a 

***  ~ — L Act 

c ly  ’’’a  2p  + C^. 
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Figure  2-2.  The  A-Matrix  for  the  Aerodynamic  Equations  of 
Longitudinal  Flight  in  State  Space  Form. 
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Figure  2-3.  The  B-Matrix  for  the  Aerodynamic  Equations  of 
Longitudinal  Flight  in  State  Space  Form. 
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In  the  longitudinal  flight  configuration  for  an  aircraft,  there 
are  two  distinct  characteristic  response  modes.  The  first  is  a short 
period  mode  which  usually  has  a natural  frequency  of  approximately  3 
radians/sec  in  comparison  to  the  second,  called  the  phugoid  mode,  which 
has  a rather  long  period  (usually  between  50  and  150  seconds).  These 
response  modes  are,  of  course,  related  to  the  eigenvalues  of  the  system 
which  in  turn  are  dependent  on  the  aerodynamic  coefficients  of  the 
aircraft. 

The  eigenvalues  are  solutions  to  the  characteristic  equation 
which  is  given  by 


C.E.  = jxl-Al  = 0 


(2-18) 


where 

x4  eigenvalue 

l4  Identity  matrix 

A = A matrix  in  Equation  (2-17) 

The  characteristic  equation  is  an  n^*^  order  polynomial  where  n is  the 
dimension  of  the  A matrix. 

For  the  fifth  order  aircraft  model  there  will  be  five  eigenvalues, 
which  may  be  real  or  complex.  There  are  three  possibilities  that  must 
be  examined. 

1.  All  eigenvalues  are  real. 

2.  Three  eigenvalues  are  real  and  two  are  complex  conjugates. 

One  eigenvalue  is  real  and  four  are  complex  (two  pairs  of 
complex  conjugates). 


3. 
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The  third  case  is  by  far  the  nx)st  coimon  with  aircraft  although  2 and  1 
do  occur.  If  the  eigenvalues  for  the  system  fit  the  third  case,  the 
characteristic  equation  can  be  written 

(A  + c^)(x^  + b2X  + C2)(a^  + b3X  + C3)=0  (2-19) 

The  second  order  factors  of  the  characteristic  equation  can  be 
represented  by 

x2  + 2cw  X + w 2 (2-20) 

n n 

where 

C = damping  ratio 

0)^  = undamped  natural  frequency 

One  second  order  factor  determines  the  response  of  the  short  period 
mode  and  the  other  provides  the  phugoid  mode  response. 

An  inspection  of  the  generalized  A matrix  (Figure  2-2)  reveals 
that  the  elements  vary  under  different  flight  conditions,  e.g.  velocity 
and  altitude.  Simply  knowing  the  aerodynamic  coefficients  for  a given 
aircraft  under  given  flight  conditions  and  being  able  to  plug  in  the 
numbers  to  compute  the  elements  of  the  A matrix  to  calculate  the 
eigenvalues  for  one  situation  leaves  much  to  be  desired  with  respect 
to  the  general  aircraft  control  problem.  In  the  design  of  the  com- 
pensator it  is  important  to  know  which  elements  of  the  A matrix 
dominate  the  response  modes  of  the  system. 

In  order  to  analyze  the  aircraft  longitudinal  response,  schemes 
based  on  empirical  results  have  been  devised  to  approximate  the  short 
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period  and  phugoid  modes.  Since  the  z derivatives  are  small,  it  is 
almost  universal  practice  to  neglect  the  aZ  state.  Experimental 
results  reveal  the  short  period  oscillation  occurs  at  virtually  con- 
stant speed  (therefore  aV=0  in  the  short  period  mode)  while  the  phugoid 
takes  place  at  virtually  constant  angle  of  attack  (therefore  Aa=0  in 
the  phugoid  mode)  [6]. 

Short  Period  Approximation 

The  short  period  approximation  neglects  the  aV  and  aZ  states, 
reducing  the  model  to  a third  order  system.  Equation  (2-21)  gives  the 
reduced  homogeneous  equations  in  state  variable  form. 


+ 

CDe) 

2u 

1 

0 

'Da' 

2u  + 

Cl- 

‘-a 

0 

+ Cl  . 
■-a 

Aa 

D0 

= 

0 

0 

1 

e 

_Dq 

1 

C 

Cm. (Cl  + Cp  ) 
a a 

0 

1 

C + 

Cm.(2p  - Cu) 
a M 

q 

X 

2u  + Cl. 

a J 

h 

mq 

2y  + Cl  . 

'-a 

(2-21) 

Clearly,  the  reduced  matrix  is  singular  and  has  rank  = 2,  therefore  one 
of  its  eigenvalues  is  equal  to  zero  and  of  little  interest.  From  this 
approximation  it  is  clear  that  the  values  of  the  6^22,  ^2H*  ®42’  ®44 

elements  of  the  A matrix  determine  the  short  period  response.  There- 
fore the  equations  which  determine  the  short  period  response  are 
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Da 


2\i  + Cl. 


- ^Lq 

2\i  + Cl. 

a 


J 

1 

^ _ Cn^(CL^  + COg) 

1 

-sf 

ly 

''I'l 

“ 2vi  + Cl. 

a ^ 

h 

(2y  + 

Cl.) 

a -* 

|Aa 

q 


(2-22) 


The  Phugoid  Approximation 

The  phugoid  approximation  is  not  as  clear-cut  as  the  short  period 
approximation.  Unfortunately,  one  cannot  break  the  A matrix  down  into 
a submatrix  that  will  give  obvious  results  as  with  the  short  period 
approximation.  Therefore  it  is  not  possible  to  calculate  the  phugoid 
component  of  the  characteristic  equation  directly  from  any  given 
elements  of  the  A matrix.  However,  the  phugoid  can  be  predicted 
fairly  accurately  from  the  aerodynamic  coefficients  of  the  aircraft. 

Recalling  that  in  the  phugoid  mode  the  angle  of  attack,  a,  is 
virtually  a constant,  hence,  Aa=0.  This  implies  zero  pitching  moment, 
so  that  pitch  equilibrium  is  always  maintained.  This  suggests  that 
the  pitching  moment  equation  should  be  dropped.  The  reduced  equations 
are  given  by  Equation  (2-23). 


DV 

0 

D0 


2w 


^Ly  * ~ ^Lq 

+ Cl_  2p  + Cl^ 


0 


AV 

q 

AO 


(2-23) 
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Since  usually, 


2y»C| 


2u>>C. 


‘^T  ""S 

V V 


Then 


^L.’  Cp  will  be  neglected, 

a V q 


The  second  of  the  three  equations  is  an  algebraic  relation,  i.e.  with 
the  preceding  approximations 


- ^ AV  + q = 0 


(2-24) 


Using  the  approximations  and  Equation  (2-24)  to  eliminate  q.  Equation 
(2-23)  can  be  shown  to  be 


DV 

De 


‘V 

2y 


(2-25) 


Recalling  that 
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DV  = ^ 

2Vg  dt 

D9  = 

2Vg  dt 

Equation  (2-25)  can  be  written 


V Ct 

r *1 

e Wg 

av 

pc 

pF 

AV 

AG 

Ae 

0 

PC 

From  Equation  (2-26)  u)_  and  r,  can  be  determined. 

''n  " 


(2-26) 


^ Ve^e 

yc 


C 


P 


_L  ^ 

2/2 


(2-27) 

(2-28) 


With  the  knowledge  of  how  ; and  u)_  vary  with  the  aerodynamic 

^ P 

coefficients,  the  elements  of  the  A matrix  which  are  dominantly  depen- 
dent on  these  coefficients  can  be  determined  by  inspection.  An  exam- 
ination of  the  A matrix  reveals  9^3  and  a2i  are  dominated  by  C„  and  p; 
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therefore  they  effect  u)„  . The  a,,  element  is  dominated  by  Ct  and 

P y 

ul  therefore  it  effects  Cp. 

Knowledge  of  which  elemeats  of  the  A matrix  effect  the  short 
period  and  phugoid  response  modes  enables  one  to  construct  a model 
which  will  respond  to  some  desired  performance  criteria.  Because  of 
the  indefiniteness  of  the  dominate  phugoid  elements,  a trial  and  error 
technique  will  be  required  to  get  the  exact  desired  phugoid  response. 

Perfomance  Criteria 

Much  research  has  been  done  in  the  area  of  aircraft  handling 
qualities  to  ascertain  which  vehicle  parameters  influence  pilot 
acceptance.  In  investigating  the  handling  qualities  related  to  longi- 
tudinal dynamics  the  problem  is  usually  separated  into  two  parts, 
associated  with  the  short  period  response  and  the  phugoid  response. 
Attempts  are  then  made  to  correlate  pilot  opinion  with  various  para- 
meters. 

First  consider  the  phugoid  response.  For  conventional  fixed  wing 
airplanes  the  phugoid  period  is  very  long  and  not  a significant  factor 
in  pilot  rating,  however  the  phugoid  damping  is  important  [6].  As  the 
damping  decreases,  the  pilot  must  devote  more  attention  to  controlling 
the  low  frequency  motion.  One  study  by  F.  O'Hara  using  simulators 
showed  that  a damping  ratio  of  .2  was  rated  acceptable  and  one  of  .6 
was  rated  good  by  the  pilots  participating  in  the  study  [6]. 

The  parameters  that  effect  the  pilot  ratings  of  the  short  period 
response  are  more  complex  than  the  phugoid  response.  Substantial 
disagreement  among  results  based  on  simply  correlating  pilot  ratings 
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with  short  period  damping  and  natural  frequency  has  resulted  in  a search 
for  more  meaningful  parameters.  One  such  parameter  was  derived  by 
noting  that  the  pilot's  opinion  of  an  aircraft's  handling  qualities  is 
influenced  by  the  vehicle's  response  to  control  inputs.  This  depends 
on  both  the  poles  and  zeroes  of  the  system  transfer  function.  An 
important  transfer  function  is  the  approximate  one  relating  pitch  rate 
response  to  elevator  angle  input  [6].  This  transfer  function  in  the  s 
domain  is  given  by  Equation  (2-29) 


A6g(s)  lyCs^  + 2C5pU)p^^S  + “n^p 


(2-29) 


To  obtain  the  impulse  response  of  Equation  (2-29)  let  A6g(s)  = 1 and 
take  the  inverse  Laplace  transform  as  given  by  Equation  (2-30) 


q(t)  = ^ 


2ct 


+ 1 


“n  “n 


sin(a)t-(j)) 


where 


and 


(2-30) 


Equation  (2-30)  shows  that  the  phase  and  magnitude  of  the  response  is 
determined  by  q and  L^/mVu)^.  These  coefficients  have  been  identified 
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a.  import  t parameters  in  the  longitudinal  handling  qualities 
criteria  [6]. 

Shomber  and  Gertsen  noted  the  importance  of  L^mv  ir.  -leir  study 

£5].  They  argued  that  when  the  aircraft's  normal  acceleration  change 

per  unit  angle  of  attack  (N^  = is  less  than  15  g/rad,  the  pilot  is 

concerned  with  controlling  the  aircraft's  flight  path.  The  magnitude 

of  the  flight  path  curvature  due  to  elevator  deflection  is  approximately 

(L^/mv)Aa.  However  when  > 15  g/rad,  the  stress  experienced  by  the 

pilot  due  to  normal  acceleration  forces  causes  the  pilot  to  be  more 

concerned  with  controlling  the  normal  acceleration  than  controlling  the 

flight  path.  Therefore  Shomber  & Gertsen  derived  two  sets  of  parameters. 

The  first  (for  N < 15)  is  a function  of  L /mv,  u-  , and  c-r,  and  the 
a ' a ^sp 

second  (for  > 15)  is  a function  of  , and 

Figures  2-4  and  2-5  show  iso-opinion  curves  based  on  the  use  of 
these  parameters.  The  data  used  to  form  these  curves  is  based  on  pilot 
ratings  from  experiments  using  simulators  and  variable  stability  air- 
craft. The  solid  lines  represent  curves  of  constant  pilot  rating  as 

the  values  of  L /mvaj^  or  N , as  the  case  may  be,  and  are 

sp  “ "sp 

varied.  The  regions  of  satisfactory,  acceptable,  and  unacceptable 

handling  qualities  are  indicated. 

From  Figures  2-4  and  2-5  it  is  evident  that  ideally  L /u„,  and 

a n 

C should  be  constant  for  all  flight  conditions.  Since 

L = l/2pV^SC,  = KpV^ 
a L 


a mg 


J 


and 


;t  r.l 
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then  the  compensated  system  would  have  to  meet  the  following  conditions: 
1.  c is  not  function  of  p or  V. 


For 

N 

< 15 

a 

U) 

n 

= kpV 

(2-31) 

For 

N 

> 15 

a 

“n 

= kpV^ 

(2-32) 

However,  from  a practical  viewpoint,  the  previous  three  rigid  conditions 
do  not  have  to  be  met  as  long  as  the  compensated  system  response  remains 
in  the  satisfactory  region. 

If  it  is  assumed  that  the  compensated  system  will  employ  state 
feedback  of  the  form 

X = [A  + BK]x  + Bu  (2-33) 

where  the  K matrix  is  a constant  feedback  matrix,  then  the  characteristic 
response  of  the  compensated  system  will  be  determined  by  the  matrix 
[A  + BK].  Since  the  elements  of  both  the  A and  B matrices  are  depen- 
dent on  the  aerodynamic  coefficients  for  the  system,  the  compensated 
system  response  will  be  a function  of  the  A and  B matrices.  Recalling 
that  the  short  period  response  of  the  open  loop  (uncompensated)  system 
is  dependent  mainly  on  the  Aa  and  q states,  it  will  be  assumed  that 
the  compensated  system's  short  period  response  is  also  dependent  only 
on  these  states.  Using  Equation  (2-22),  the  second  and  fourth  rows  of 
the  B matrix  (Figure  2-3)  and  the  definition  for  Da,  Dq,  and  q. 

Equations  (2-34)  and  (2-35)  are  determined. 


c + Cl- 


2w  - CLq 

2u  + Cl* 
a 


Cn,^(CL^  + CDg) 
r 

2Ve‘ 
r 

a 2w  + Cl  . 

L Lqi  J 

iyc[ 

Cin^(2p  - CLq) 


(2-34) 


c(2vi  + C|_.) 


c(2p  + Cl  .) 

a 


(2-35) 


Iw  1 

1 (2p  + Cl  .) 

a -* 

■y'  . 

where  Kc  = C„  $111(0^)+  1 
b nr  ' T' 
a 

and  , K2,  K^,  and  are  constants  that  are  functions  of  the 
elements  of  the  feedback  matrix.  Making  the  following  assumptions, 

2p»C|_. 


2ij»C, 


I..  = — 


-h. 

V = 


Cr,  , c,  C , C , C , C,  are  all  constants  then  Equations  (2-34)  and 
m m*  m„  L 
e a a q a 

(2-35)  become 


F-  r 


Although  Equations  (2-42)  and  (2-43)  do  not  quite  meet  the  con- 
ditions for  the  ideal  system,  if  the  compensator  is  designed  so  that 
at  a median  velocity  and  altitude  the  compensated  system's  response  is 
well  within  the  satisfactory  region,  then  for  a certain  limit  of  flight 
conditions  the  system  response  will  remain  in  the  satisfactory  region. 

Since  the  location  of  the  zero,  in  Equation  (2-29)  is  an 

important  parameter  in  the  performance  of  the  aircraft,  it  will  be 
useful  to  determine  which  elements  of  the  A matrix  determine  the  zero 
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location.  Using  the  short  period  approximation  matrices,  the  transfer 
function  can  be  determined  by 


= C[sI-A]"^B 


A6«(s) 


where 

C = [0  1] 
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CM 

42 

"44 

K]b22  + K2b23 

B = 

= 

K3b42  + K4b43 

(2-44) 

(2-45) 

(2-46) 


(2-47) 


From  Equations  (2-45),  (2-46),  and  (2-47),  Equation  (2-44)  can  be  shown 
to  be 


q(s)  b2'(s  - 322^  ” ^^42 

A^e^s)  s2  + 2^UnS  + 0)^2 


(2-48) 


If  b2^»bi",  then  Equation  (2-48)  reduces  to 


q(s)  b2'(s  - ^22^ 
ACe(s)  s2  + 2ca)^S+ 


(2-49) 


Equation  (2-49)  shows  the  zero  location  is  determined  by  the  a22  element. 


i 

i 


29 


The  3^2  element  can  be  expressed  as 


'22 


c 2u  + C|  •• 

^ni 


If  it  is  assumed  that 


2vi>>C, 


Cl  »Cjj 


then  Equation  (2-50)  becomes 
VCl 


ao.. 


Since, 


22  wc 


2L 


and 


y = 2m 


pSc 


(2-50) 


(2-51) 


then 


(2-52) 


Therefore  the  zero  is  located  at  s = L /mv  which  is  the  same  result  as 

a 

given  in  Equation  (2-29). 


With  the  knowledge  of  how  the  elements  of  the  A matrix  effect  the 
system  response,  a method  for  developing  a satisfactory  model  can  now 
be  formulated. 
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1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 

9. 

10. 

11. 


Through  simulation  or  other  means  choose  a median  flight 
condition  for  which  the  compensator  will  be  designed. 

Pick  a desired  c and  for  the  phugoid  and  adjust  the  a^^* 
32^  and  a-ji  elements  for  the  desired  response. 

Determine  N at  the  chosen  velocity  and  altitude  from  the 

“ L 

aerodynamic  coefficients  using  N = — . 

a mg 

If  N < 15,  then  determine  L /mv  and  calculate  an  and 

a a n^p 

a which  place  the  intersection  of  L /mvu  and  c._  in 
sp  o n^p  sp 

the  center  of  the  satisfactory  region  in  Figure  2-4. 

If  N >15,  then  determine  an  un  and  Ccn  which  will  place 

a "sp 

the  intersection  of  N and  c in  the  center  of  the 

“ sp  P 

satisfactory  region  in  Figure  2-5. 


Using  the  desired  values  of  and  , calculate  the  short 

P sp 

period  characteristic  equation  of  the  model 

C.E.  = s^  + 2c__u)„  s ^ 

sp  n^„  n._ 

^ sp  sp 


Set  ap.  = 1 

Do  not  change  a22.  Therefore  822  = 322  • This  is  done  so 

^m  p 

the  zero  location  in  Equation  (2-29)  is  not  changed. 

Determine  a^^  by  a.^  = - (a22  + 2c  ^ ) 

m m m ^ sp 

Determine  a^2  by  a. 2 = (^22  ®44  ^ ‘ “n^' 
m m m m 

Simulate  the  model  to  insure  that  the  model  has  the  desired 
response. 


With  the  knowledge  of  how  to  design  a model  that  will  perform  as 
per  the  performance  criteria,  the  task  at  hand  is  to  design  a compen- 
sator which  will  cause  the  aircraft  to  have  the  same  response  as  the 


model . 


III.  DEVELOPMENT  OF  COMPENSATOR  DESIGN  ALGORITHMS 


Once  a model  is  designed  that  will  perform  as  per  the  specified 
performance  criteria,  the  problem  at  hand  is  the  design  of  a compen- 
sator that  will  force  the  plant  to  respond  in  the  same  manner  as  the 
model.  There  are  certain  physical  restraints  with  the  aircraft  problem 
that  must  be  kept  in  mind.  Consider  the  block  diagram  for  the  mathe- 
matical aircraft  model  as  given  in  Figure  3-1.  First,  it  is  assumed 
that  nothing  within  the  dotted  outline  of  Figure  3-1  is  accessible. 

In  other  words,  the  A and  B matrices  cannot  be  altered  directly  since 
they  are  determined  by  the  configuration  of  the  aircraft  and  its 
environment.  Second,  it  is  assumed  that  the  states  of  the  aircraft 
{x  vector)  are  available  which  means  that  there  are  devices  that  can 
accurately  measure  the  velocity,  angle  of  attack,  pitch  angle,  pitch 
rate,  and  altitude.  Since  the  only  available  information  in  the 
system  are  the  inputs  and  the  states,  the  compensator  C(s)  must  be  of 
the  form  given  in  Figure  3-2. 

R.  T.  Curran  developed  an  algorithmic  approach  for  the  design  of 
the  compensator  utilizing  his  concept  of  "equicontrollabil ity"  [4]. 

Definition:  Given  the  system 

X = Ax  + Bu^  ( 3- 1 ) 

where  A has  dimension  n x n and  rank  n 

B has  dimension  n x m and  rank  m 
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Figure  3-1.  System  Block  Diagram  Representation  of  Mathematical 
Aircraft  Model. 
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the  system  is  controllable  if  and  only  if  there  is  a p_<n  such 
that 

Cp  = [B,  (3-2) 

has  full  rank  n = dimension  of  A and  furthermore  the  pair  (A,B) 
has  controllability  index  p if  p is  the  smallest  integer  that 
satisfies  Equation  (3-2). 

Definition;  The  pair  (A,B)  is  equicontrol lable  if  n = pm  where 
p is  the  controllability  index  of  (A,B),  A is  n x n, 
and  ” is  n X m. 


Note,  the  definition  of  equicontrollability  implies  that  the  first 
n = pm  columns  of  Cp  are  independent. 

Based  on  Luenberger's  work  on  canonical  forms  for  multivariable 
systems  [7],  Curran  showed  that  if  a system  is  equicontrol lable  then 
there  exists  a similarity  transformation  T such  that 

m n-m 

T'^AT  = (3-3 

[ X J m 


and 


T'^B  = 


LU 


n-m 


(3-4) 


The  restriction,  that  for  a system  to  be  equicontrol lable,  (n=pm) 


seems  to  be  severe  and  therefore  of  little  practical  usefulness.  Indeed, 
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most  systems  will  not  meet  this  structural  constraint.  To  overcome 
this  problem,  Curran  proposes  to  add  states  to  the  system  that  will 
alter  the  structure  of  the  system  but  not  it's  transfer  function.  Such 
states  will  have  to  be  controllable  but  not  observable.  Curran  proves 
the  following: 

Given  the  system  (A,B)  with  controllability  index  p,  A is  n x n, 

B is  n X m.  Then  if  n < pm,  there  exists  matrices  J and  L and  a number 
N such  that  the  N-dimensional  system 


is  equicontrollable  if  and  only  if  N = m(p  + k)  for  some  integer  k ^ 0. 

The  manner  in  which  the  additional  k states  are  connected  to  the 
original  system  will  determine  the  J matrix  and  the  pole-placement  of 
the  additional  states  will  determine  the  L matrix. 

Before  pursuing  the  determination  of  the  J and  L matrices  it  will 
prove  useful  to  examine  the  structure  of  an  equicontrollable  system  in 
canonical  form.  Consider  the  following  example. 

Example  3.1 

Given  the  equicontrollable  canonical  system  (A,B)  with  control - 


J 


lability  index  p,  A is  n x n,  B is  n x m 
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^4  ^ 

Xg  ~ ^2^—’—^ 

^6 "" 


FroTi  the  state  equation  a system  block  diagram  can  be  readily  drawn  as 
shown  in  Figure  3-3. 

There  are  3 strings  of  integrators,  each  consisting  of  two 
integrators.  In  general,  for  an  equicontrollable  system  in  canonical 
form  there  will  be  m strings  of  integrators  of  length  p.  In  other 
words,  each  of  the  inputs  controls  an  equal  number  of  states. 

Curran  shows  that  given  the  system  matrices  (A,B),  where  A is  n x n 
and  B is  n X m and  full  rank,  then  there  exists  a similarity  transfor- 
mation T such  that 


AT  = 


m n-m 

I 

0 I I 


X 


n-m 


m 


(3-5) 


B = 


0 n-m 

W m 


(3-6) 


if  and  only  if  the  last  n columns  of  the  controllability  matrix 


C = [A^'^B,  a"‘^B,...AB,B] 


are  linearly  independent. 
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f2(x,u) 


/ 


Figure  3-3.  Block  Diagram  Representation  of  State  Equations  for 
the  Equicontrollable  System  Given  in  Example  3.1. 
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Once  the  system  has  been  transformed  into  the  form  of  equati^.is  i-5) 
and  (3-6)  the  states  can  easily  be  added  so  the  system  is  in  equicon- 
trollable  form.  Example  3.2  illustrates  this  process. 

Example  3.2 

Given  the  system  matrices  (F,G)  in  the  form  of  Equations  (3-5) 
and  (3-6)  where  F is  5x5,  G is  5x3 


1 

0 

0 

0 

0 

1 

1 

0 

0 

0 

1 

0 0 0 0 1 

0 0 0 

^1 

G = 

-91- 

^2 

— 92- 

^3 

_ — 93 — 

then 

Xi  = X4 
X2  = X5 

X3  = fi(2L)  + g|(y.) 

X4  = f2(x)  + QzCm.) 

X5  = f3(x)  + g3(u) 

and  the  system  block  diagram  will  be  the  form  as  shown  in  Figure  3-4. 

From  the  block  diagram  it  is  obvious  that  the  addition  of  one  integrator 

will  make  the  system  equicontrol lable.  The  state  (xg)  should  be  added 
such  that 

^6  ~ ^3  ^ aXg 


(3-7) 


40 


where  a is  the  pole  location  of  the  added  state.  The  augmented  system 
matrices  (F,6)  are  now  in  the  form 


1 

o 

1 

iO 

o 

u. 

and  G = 

G 

o 

1 

1 

1 

1 

o 

» 

t 

1 

O 

o 

(3-8) 


Hence, 


J = [0  0 1 0 0]  (3-9) 

and 

L = [o]  (3-10) 


Once  the  system  is  expanded  such  that  it  is  equi control! able  then 
the  system  can  once  again  be  transformed  so  the  final  system  matrices 
are  in  equicontrol lable  canonical  form. 

So  far  it  has  only  been  shown  that  there  exists  a similarity 
transformation  matrix  T that  will  transform  the  system  matrices  into  the 
desired  form.  The  process  for  generating  this  transformation  matrix  T 
is  straightforward  and  is  outlined  in  the  following  steps. 

1.  Given  the  system  matrices  (A,B)  with  controllability  index  p, 
where  A is  nxn  and  B is  nxm,  generate  the  matrix  P from  the 
last  n columns*  of  the  partial  controllability  matrix  C where 

*(the  last  n columns  must  be  taken  so  the  transformed  system  will 
be  of  the  form  given  in  Equation  (3-5).) 
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(3-11) 


2.  Compute  the  matrix  Q such  that 


Q = P 


-1 


(3-12) 


3.  Take  the  first  m rows  of  Q to  form  the  mxn  matrix  E and  par- 
tition it  such  that 


E = 


(3-13) 


m-r 


where  r = n - (p-1 )m 
4.  Form  the  matrix  S such  that 


E 

EA 


eaP- 


e^aP-^ 


(3-14) 


5.  Determine  the  transformation  matrix  T where 


T = S 


-1 


(3-15) 


1 


Note  that  the  above  algorithm  works  for  all  systems  that  meet  the  require- 
ment that  the  last  n columns  of  the  controllability  matrix  are  linearly 
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independent.  For  an  equicontrollable  system  n = pm  and  m-r  in  Equation 
(3-13)  would  equal  zero. 

The  purpose  of  augmenting  .the  system  to  achieve  equicontrollability 
and  then  transforming  it  to  a canonical  form  is  that  the  structure  of 
the  transformed  system  can  be  utilized  to  easily  find  a state  feedback 
matrix  K that  will  alter  the  pole  and  zero  locations  of  the  original 
system.  The  problem  is  outlined  below. 

Given  the  plant  and  model  state  equations 


^ + B u 

-P  P-P  P- 


jc  = Ax  + B u 
-m  m-m  rrr- 


(3-16) 

(3-17) 


where 


dimension  of  A and  A = nxn 
p m 

dimension  of  and  B = nxm 
p m 

determine  the  feedback  matrix  K such  that  the  transfer  function  of  the 
compensated  plant  equals  the  transfer  function  of  the  model  as  given  by 
Equation  (3-18). 

[SI  - (Ap  * BpK)r'8p  = [SI  - A^r'B^  (3-18) 

The  process  for  determining  the  matrix  K involves  parallel 
operations  on  the  plant  and  model  state  equations  beginning  with  the 
original  systems  given  by  Equations  (3-16)  and  (3-17).  An  overview  of 
the  algorithm  is  given  on  the  following  page. 


w' 


i 

! 
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1.  Transform  the  plant  and  model  to  a canonical  form  as  given  by 
Equations  (3-5)  and  (3-6). 

2.  Add  states  to  the  plant  arili-.niodel  as  outlined  in  Example  3-2. 

3.  Transform  the  augmented  plant  and  model  systems  to  equicon- 
trollable  form. 

4.  If  the  plant  and  model  similarity  transformations  used  in 
1 and  3 are  equal,  then  solve  directly  for  the  feedback 
matrix  K. 

5.  Transform  the  compensated  plant  back  to  the  original 

coordinates.  ^ 

Note  from  4,  the  feedback  matrix  K can  be  solved  directly  if  the 
plant  and  model  similarity  transformations  are  equal.  It  will  be  shown 
later  that  for  aircraft  longitudinal  equations  the  similarity  trans- 
formations will  always  be  the  same  if  care  is  taken  in  constructing  the 
model.  For  systems  that  do  not  meet  this  requirement,  Curran's  algorithm 
completes  the  procedure  for  developing  the  compensator. 


The  Algorithm  in  Detail 

Given  the  state  equations  for  the  plant 


= Ax  + 
-P  P-P 


(3-19) 


transform  the  system  to  a canonical  form  so  the  additional  states  can 
be  added.  Let 


(3-20) 


Substituting  Equation  (3-20)  into  Equation  (3-19)  gives 
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n = A X + B U 

PpP  P P^“p  P— 

(3-21) 

(3-22) 

X 

= Ax  + B u 

-p  p-p  p- 

(3-23) 

where 


A 


P 


-1 


m (n-m) 


(n-m) 


m 


(3-24) 


and 


e 


n-m 


(3-25) 


Augment  the  above  equations  such  that  the  new  system  is  equicontroHable 
(N=pm)  as  shown  in  Figure  3-5  and  given  by  Equation  (3-26).  Let 


(3-26) 


then  the  state  equations  become 


3-5. 
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I- 


(3-27) 


= FJL  + G u 

-p  p-p  p- 


(3-28) 


where 


(n  + r)  = N = pm 


The  J and  L matrices  are  determined  by  the  method  given  in  Example  3.2. 

Next  transform  the  augmented  system  to  the  equi control! able 

canonical  form  via  similarity  transformation  T . Now  the  system 

^2 

equations  are 


where 


— * -p-  * — * . K * 

x„  = F„  x„  + u 
-P  P-P  P - 


(3-29) 


_ _ _ * 
Xd  = x„ 

-P  p^-p 


(3-30) 


— * -1— 

F = T FT 
P P2  P P2 


(3-31) 


G = T G 
P P2  P 


(p-l)m 

I 1 (p-l)m 


(3-32) 


(3-33) 
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r * 


Cp-1  )ni 


(3-34) 


Once  the  plant  has  been  augmented  and  transformed  to  equicontrol- 
lable  canonical  form,  repeat  the  same  process  with  the  model.  Once  the 
system  model  has  been  augmented  and  transformed  it  will  be  of  the  form 


m -m 


^m*M 


m (p-l)m 

6 I I (p-l)m 

m 


e 


Cp-l)m 


I 


m 


(3-35) 


(3-36) 


(3-37) 


If  T_  = T„  and  = T„  , then  the  K matrix  can  be  determined  directly 
P]  ni]  ¥2  '^2 

as  shown  in  the  following  paragraph. 

Given  the  equi controllable  canonical  system 


let 


— ★ _ ★_  ★ _ * 


it 

u = u + K 2^ 
- - -P 


Substituting  Equation  (3-39)  into  Equation  (3-38)  gives 


^ (^p  Gp  K ))^  + Gp  U. 


(3-38) 

(3-39) 


(3-40) 
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The  goal  up  to  now  has  been  to  make  the  ^'•ansfer  function  of  the 
compensated  plant  equal  to  the  transfer  function  of  model  such  that 

-1  . -1 

[si  - (A  + B K)]  B = [si  - A ] B (3-41) 

^pp'-^p*-  m-'m  ^ ' 

Utilizing  the  fact  that  the  similarity  transformations  used  for  the 
plant  and  model  are  the  same  and  Gp  = (it  will  be  shown  that  this 

implies  B|^=Bp),  the  problem  is  now  to  find  a K*  matrix  such  that 

[si  - (Fp*  + GpT)]-"'  = [si  - (3-42) 

resulting  in 


Because  of  the  structure  of  G^ 


(3-43) 


(3-44) 


where  Xp  and  are  the  mxpm  submatrices  given  in  Equations  (3-33)  and 
(3-36)  respectively.  Once  1^*  is  determined  then  the  system  is  returned 
to  the  original  coordinates.  The  process  for  returning  the  original 
coordinates  is  tedious  but  straightforward.  Since  the  plant  and  model 
transformation  matrices  were  equal,  for  convenience  the  subscripts  will 
be  dropped.  In  order  to  separate  the  added  states  from  the  states  of 
the  original  system,  T2"^  is  partitioned  such  that 


n 


r 


(3-45) 
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and 

H*  " '''21’^^  ^2?’’- 

The  state  equations  for  the-  ompensator  in  the  original  coordinates  are 

+ JT^'^x  (3-47) 

u = u + 1^121  + i^T22‘V  (3-48) 

The  block  diagram  for  the  compensator  is  given  in  Figure  3-6  and  the 
compensated  system  is  shown  in  Figure  3-2.  The  following  example 
demonstrates  the  algorithm  in  detail. 


Example  3.3 

Given  the  plant  and  model  matrices. 


-1 

2 

0 

0 

0 

Ap  = 

2 

3 

1 

1 

1 

-1 

2 

-1 

0 

2 

'-1 

2 

o' 

"0 

o' 

-1 

-2 

1 

B = 
m 

1 

1 

1 

-1 

-2 

0 

2 

The  controllability  index  p = 2.  This  means  that  N = pm  = 4; 
therefore  one  state  will  have  to  be  added  so  that  the  system  is  equi- 
control lable. 

First  the  system  will  be  transformed  to  a canonical  form  using  a 
similarity  transformation  matrix  Tp  in  order  to  determine  how  to 


Figure  3-6.  Block  Diagram  Representation  of  the  Compensator. 
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connect  the  additional  state.  Transformation  matrix  T is 

P] 

using  Equations  (3-11)  through  (3-15).  From  Equation  (3-11) 


The  matrix  P is  taken  as  the  last  n rows  of  Cp. 


P = 


2 

5 

0 


0 

1 

0 


0 

1 

2 


From  Equation  (3-12) 


.5 

0 

0 

-2.5 

1 

-.5 

0 

0 

.5 

From  Equation  (3-13) 


E = 


0 

1 


and 


generated 


From  Equation  (3-14) 


.5  0 


-2.5  1 

-.5  1 


From  Equation  (3-15) 


= S- 


2 0 0 
1 0 1 


■8  -2  2 


Using  Tp  the  original  plant  is  transformed  to  the  new  coordinate 
system  ^ where 


-P  " ^P^-P 


Using  Equation  (3-24) 


Ap  = Tp,  'ApTp,  .-5-3  0 

-3  -2  4 


0 0 1 ] 


(3-49) 


Using  Equation  (3-25) 


0 0 

B.  = = 1 0 


'P  'Pi  “P 


1 1 


(3-50) 


Using  the  process  as  demonstrated  in  Example  3.2  and  the  state  equations 
from  Equations  (3-49)  and  (3-50)  the  block  diagram  for  the  system  can 


'I 
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be  drawn  as  shown  in  Figure  3-7.  From  the  block  diagram  it  is  obvious 
that  the  added  state  will  have  to  satisfy 

z-j  = X2  (3-51) 

The  pole  of  this  added  state  is  placed  at  s = -10.  Therefore 

J = [0  1 0]  (3-52) 


and 

L = [-10] 


The  augmented  system  is  now 

Xp " ^pi^p 


(3-53) 


(3-54) 


- 

I 

■*  — 

^1 

0 0 1 1 0 

^1 

1 

o 

o 

1 

X 

y 

-5  -3  0 ' 0 

1 

^2 

+ 

1 0 

h 

-3  -2  4 ' 0 

--  -1-  - 

^3 

1 1 

0 1 0 '-10 

1 

.^1. 

_0  0_ 

(3-55) 


Next,  the  equi controllable  system  is  transformed  to  the  equi con- 
trollable canonical  form  via  similarity  transformation  matrix  T . 

P2 

Matrix  T is  generated  by  the  same  algorithm  as  matrix  T . 

P2  P] 


110  0 
10  0 1 0 
0 0 11 
10  0 0 


(3-56) 


F„*  = T "^F 
P P2  P P2 


0 1 0 

0 0 1 

-5-13  0 

2 15  4 


(3-57) 


r ■''g  = 
P2  P 


0 0 
0 0 
1 0 
0 1 


(3-58) 


The  same  process  is  now  performed  on  the  model.  Since  the  solution  is 

dependent  on  the  necessity  that  T " " ^o  perform 

•^11  1 ^1 


the  transformation. 


2 0 0 

1 0 1 

•8  -2  2 


(3-59) 


If  the  transformed  model  is  of  the  desired  form,  the  algorithm  is 
satisfied.  If  the  transformed  model  is  not  of  the  desired  form,  then 
there  is  no  solution  and  the  algorithm  should  be  terminated.  The  model 
is  then  transformed  using  T ^T 

mi  p^ 


A = T "^A  T 
m m^  m m^ 


0 

0 

1 

-20.5 

-4 

-2.5 

-12 

-2 

-1 

0 0 

-,  B = T -■'b  = 1 0 

m m^  m 

Ll  1. 

(3-60) 


1 
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The  model  is  now  augmented  so  the  system  is  equicontrollable.  It  is 
important  that  the  pole  location  of  the  added  state  for  the  model  is 
the  same  as  that  of  the  plant  t.o  insure  that  the  similarity  transfor- 


The  augmented  system  is  now 


X 

-m 


F X 


tirm 


+ G u 
m— 


(3-61) 


(3-62) 


T is  determined  and  the  augmented  model  is  transformed  to  equicon- 
^2 

troll  able  canonical  form. 


(3-63) 


(3-64) 
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’0  0“ 

0 0 

1 0 

(3-65) 

_0  1_ 

Since  Tp^ 

^m] 

and  Tp2 

= feedback  matrix 

can  be  determined 

directly. 

From 

Equation  (3-57) 

■-35  - 

5 -13 

on 

Xp  = 

(3-66) 

P 

. 12 

2 15 

4J 

From  Equation  (3-64) 

X = 
m 

■-60.5 

-20.5 

-16.5  -2.5] 

(3-67) 

. 28.5 

8.5 

13.5  I.5J 

Using  Equation  (3-44) 


x = - Xp  = 


[■-25.5  -15.5  -3.5  -2.5‘] 

L16.5  6.5  -1.5  -2.5J 


(3-68) 


To  transform  the  compensator  back  to  the  original  coordinates  the 
matrices  L,  ^^21~^^l'^*  ^^22~^  determined.  First 

is  partitioned  as  per  Equation  (3-45)  to  determine  T^^ and  ^22”^- 


~0 

0 

0 

1 1 

1 

0 

0 

0 

1 

0 

-10 

_0 

-1 

1 

10_ 

(3-69) 
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From  Equation  (3-69) 


(3-70) 


Performing  the  required  matrix  operations  the  results  are: 


L = [-10J 

JT^'^  = [-2.5  1 -.5] 


The  compensator  in  matrix  form  is 

= [-10]z^  + [-2.5  1 -.5])^ 


(3-71) 

(3-72) 

(3-73) 

(3-74) 

(3-75) 


(3-76) 
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u = 


1 0 1-4 

I 

0 1 I 1 


-3.5 

-1.5 


.5 

-.5. 


'0' 

x„ 

_o. 

■1 


(3-77) 


Note  from  Equation  (3-77)  that  the  output  of  the  compensator  ^ is  not 
a function  of  the  additional  state  Zp  but  only  of  the  input  and  the 
system  states.  Therefore  there  is  no  need  to  generate  and  the 
system  can  be  compensated  by  pure  state  feedback.  The  final  form  of 
the  compensated  system  is  shown  in  Figure  3-8  and  the  state  equations 
for  the  compensated  system  are 


where 


Since 


and 


then 


i^p  = V^^p  + ^py 


(3-78) 


Ap  + BpK  = 


-12  0 
-1  -2  1 

1 -1  -2 


= A 


m 


A + B K = A 

ppm 


Bn  = B„ 
p m 


Isl  - (A„  + B^K)j"^B„  = [si  - Aj'^B„ 
p p -'p*-  m m 


(3-79) 
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Example  3.3  was  a contrived  problem,  however  it  illustrates  a 
fundamental  flaw  in  Curran's  algorithm.  If  T had  been  generated  using 

mi 

the  method  given  by  Equations  C3-11)  through  (3-15),  T_  would  not  have 

1 

equaled  T which  would  have  created  considerably  more  work  in  solving  the 

n 

problem.  An  alternate  method  which  will  preclude  the  possibility  of 
this  happening  is  proposed  below. 


Proposition  1 

Consider  the  system  pair  (A,B)  such  that 


X = Ax  + Bu 


where  A has  dimension  n x n and  B has  dimension  n x m and  both  (A,B) 
are  full  rank  and  the  system  is  controllable  with  controllability  index 
p.  If  n ^ pm  then  add  r states  to  the  system  such  that  N = n + r = pm 
and  the  augmented  system  is  controllable  with  controllability  index  p. 

Proposition  1 removes  the  need  for  transformation  matrix  T-].  However, 

there  is  still  the  possibility  that  Tq  will  not  equal  T . The 

^2  '"2 

designer  of  course  could  simply  let  = Tp^  as  was  done  in  Example  3.3 

with  T-,  . However,  it  would  be  nice  if  there  were  some  criteria  for 
1 

insuring  that  = Tp^.  An  examination  of  the  similarity  transformation 
and  the  expected  results  will  be  useful  for  this  purpose. 

Curran  showed  that  given  the  equi controllable  system 


X = Ax  + Bu 
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i 


tit^. . .t. 


0 • ••  0 * “0 
• • 

• • 

Q...  0-  .■] 


(3-83) 


therefore  from  (3-83),  the  last  m columns  of  T must  equal  the  columns 


of  B. 


T = 

1 

T 

r ■ 

'n-m* 

’m 

I 

(3-84) 
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Using  the  results  from  Equation  (3-84)  and  Equation  (3-80) 

TT"^AT  = TA*  (3-85) 


^0--0  1 1 0 0" 

0---6  1 • ■ : 

• 1 ■ 

c 

2 • ■ • ^ 

n 

t 

^...1 

n-m  ^ 

>1 . . .t 

m 

= 

t 

^ • t 

n m t 

n-m 

*1 . . . t 

m 

1 — 

1 • • 

1 

1 

1 — 

1 °- 
1 • 

1 

1 “* 

_ 

( — 

X 

L_ 

(3-86) 


Equation  (3-84)  shows  that  for  the  plant  and  model  transformations  to 

be  equal,  the  input  matrices  of  the  two  systems  must  be  equal.  Hence 

a necessary  but  not  sufficient  criterion  for  T = T is 

P2  '‘^2 

Bm  = Bp  (3-87) 

Equation  (3-86)  shows  that  the  remaining  terms  of  the  similarity  trans- 
formation depend  on  the  elements  of  the  A and  B matrices. 

If  the  system  is  of  a particular  form.  Equation  (3-86)  will  be 
useful  in  showing  which  elements  of  the  A matrix  determine  the  trans- 
formation matrix. 

THEOREM  I 

Given  the  equi controllable  system  (A,B)  with  controllability 
index  p = 2 where  A is  dimension  n x n and  B is  dimension  n x m and 
full  rarik. 

If  B has  n - m zero  rows  then  the  similarity  transformation 
matrix  T that  transforms  the  system  to  equicontrollable  canonical  form 
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will  depend  only  on  the  B matrix  and  the  row  elements  of  the  homogeneo.^ 
states  in  the  A matrix. 


Proof 


Gi  ven 


A = 


m 

An 


m 

A. 


and  B = 


m 

e 


m 


(3-88) 


then  as  shown  in  Equation  (3-84) 
m m 


T = 


'1 


e 


m 


m 


(3-89) 


Since  T”^AT  = A*  where  A*  is  the  A matrix  in  equi controllable  canonical 
form,  then  AT  = TA*  as  given  in  Equation  (3-90). 


A-]  1 1^2 


Ti  I 6 


X ' 
^ I 


From  Equation  (3-90),  Equations  (3-91)  and  (3-92)  are  obtained 


^1^1  ^2^2  " ® 


^2^1  ""  ^1 


"2 
(3-90) 


(3-91) 


(3-92) 
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Writing  Equations  (3-91)  and  (3-92)  in  matrix  form, 
m m 

1 r 1 

m 

[D] 


where 


- - 

— 

h 

m 

h 

m 

1 

CO 

CNJ 

c 

J 

(3-93) 


m 


m m 


[D] 


In  order  for  a unique  solution  to  exist  for  T-j  and  T2  in  Equation  (3-93), 
the  matrix  D must  have  an  inverse.  By  inspection  the  matrix  D will 
have  an  inverse  if  and  only  if  A2  has  rank  = m.  Since  the  controllability 
matrix  C where 


m m 


^2^1  I ^ 
1 

A4B, ] B, 


m 


m 


has  rank  = n and  has  rank  = m then  must  have  rank  = rn.  Since 
will  always  have  rank  = m for  the  system  to  be  equi controllable  with 
controllability  index  p = 2,  then  a solution  to  Equation  (3-93)  will 
always  exist.  Therefore  the  transformation  matrix  depends  only  on  the 
B matrix  and  the  row  elements  of  the  homogeneous  states.  The  signifi- 
cance of  Theorem  I is  that  the  designer  now  ca"  change  m rows  (of  the 
non-homogeneous  states)  in  the  system's  A matrix  without  effecting  the 
similarity  transformation  matrix.  For  a system  with  controllability 
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index  p = 2 this  means  t‘"  designer  has  control  of  half  the  elements 
of  his  system  matrix  A t he  use  of  state  feedback.  Example  3.4 
utilizes  Proposition  1 a;  a Theorem  I to  solve  the  problem  presented  in 
Example  3.3. 


Example  3.4 

Given  the  plant  and  model  matrices 


“-1 

2 

o” 

~0 

0 

2 

3 

1 

1 

1 

_-l 

2 

-1_ 

_0 

2 

~-l  2 o' 

"o  0 

A = 

-1  -2  1 

1 1 

m 

_ 1 -1  -2_ 

_0  2_ 

First,  from  Proposition  1,  one  state  is  added  to  the  system  such  that 
the  augmented  system  is  controllable  with  controllability  index  p = 2. 


(3-94) 


(3-95) 
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0 0 
1 1 


0 2 
0 0 


(3-100) 


/ 

/ 


Since  the  rows  of  the  homogeneous  states  o/^the  plant  and  model  are 
equal,  then  by  Theorem  I the  transformation  matrix  for  the  model  will  be 
the  same  as  the  one  for  the  plant.  'The  transformed  model  is 


0 0 10 

0 0 0 1 

/ 

-4.5  35.5  -2.5  11.5 


.5  -19.5  -.5  -12.5 


(3-101) 


0 0 


0 0 
1 0 


0 1 


(3-102) 


The  K is  determined  in  the  same  manner  as  in  Example  3.3 


K* 


-3.5 

-1.5 


(3-103) 


Transforming  the  system  back  to  the  original  coordinates  gives  the 
same  solution  as  in  Example  3.3. 


All  of  the  algorithms  ^or  transforming  the  system  equations  to 
new  coordinates  have  imposec  le  restriction  that  the  input  matrix  B 
have  full  rank.  If  the  number  of  inputs  m is  less  than  the  order  of 
the  system,  then  this  means  the  input  matrix  B is  made  up  of  m indepen- 
dent column  vectors.  Many  systems,  including  the  aircraft  dynamics 
system,  do  not  meet  this  restriction.  Fortunately,  the  process  can  be 


extended  to  systems  that  do  not  have  input  matrices  of  full  rank. 
Proposition  2 extends  the  algorithm. 


Proposition  2 

Given  ':he  system  matrices  (A,B)  where  A has  dimension  n x n and 
B has  dime?  )n  n x r and  rank  m where  m < r then  the  matrix  B can  be 
partitioned  such  that 


where  the  submatrix  B'  has  rank  m.  The  input  vector  £ will  be 
correspondingly  partitioned 

u = 

such  that 


Bu_  = B'^'  + Rx 
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And  if  the  system  matrices  (A,B')  meet  the  criteria  for  equicontrol- 

1 ability  then  a feedback  matrix  K can  be  determined  such  that 

[SI  - (Ap  + Bp-K)]-'Bp  = [si  - (3-106) 

To  summarize  what  has  been  developed  so  far, 

1.  A system  is  or  can  be  made  equi control! able  if  the  last  n 
columns  of  its  controllability  matrix  are  independent. 

2.  If  the  similarity  transformations  that  transform  two  systems 
to  equi control  1 able  canonical  form  are  equal  then  a feedback 
matrix  can  be  determined  directly  so  that  one  system 
transfer  function  can  be  made  to  equal  the  other. 

3.  A necessary  but  not  sufficient  condition  for  the  similarity 
transformations  to  be  equal  is  that  the  systems’  input 
matrices  must  be  equal. 

4.  If  a system  is  equicontrollable  with  a controllability  index 
p=2  and  its  input  matrix  has  n-m  zero  rows,  then  the  simi- 
larity transformation  matrix  that  will  transform  it  to 
equicontrollable  canonical  form  is  dependent  only  on  the 
input  matrix  and  the  row  elements  of  the  homogeneous  states. 

An  Alternate  Method 

As  stated  previously,  the  problem  is  to  find  a matrix  K such  that 
[si  - (Ap  * BpK)r'Bp  = [si  - A„]-'b„ 

If  it  is  assumed  that  B = B , then  the  problem  reduces  to  finding  a K 

pm 


matrix  such  that 


Then, 

BpK  = tAn,  - Ap]  (3-108) 


Multiplying  both  sides  of  Equation  (3-108)  by  BpT 

IbJbJK  = bJiA^  - AJ  (3-109) 

*■  p p p m p 

If  Bp  has  dimension  n x m and  rank  m then  [Bp^Bp]  has  rank  m and 
dimension  m x m,  therefore  [Bp^Bp]  has  an  inverse, 

K " - Ap]  (3-110) 

On  the  surface  Equation  (3-110)  looks  fantastic  since  the  only  ■! 

I 

J 

obvious  restriction  for  determining  the  K matrix  is  that  the  matrix  Bp 
have  full  rank.  Unfortunately,  however,  there  are  more  obscure  restric- 
tions. First  from  Equation  (3-108),  the  rank  of  BpK  is  less  than  or 
equal  to  the  rank  of  Bp,  therefore  the  rank  of  [Aj^^-Ap]  must  be  less  than 
or  equal  to  the  rank  of  Bp.  There  are  other  structural  restraints  on  A^^ 
which  are  dependent  on  Bp.  Substituting  Equation  (3-110)  into  Equation 
(3-107)  gives  the  following  necessary  condition  for  a solution  to  the 
problem.  For  a particular  model  A^^,  there  is  a feedback  matrix  K such 
that 

[A„  + B„K]  = A 
■■  p p m 

if  and  only  if 

Ap  ^p(^p  ^pJ  ^p^[A|^  ■ Ap]  = A^  (3-111) 


1 

J 


Example  3.5  reworks  Example  3.3  using  Equation  (3-110). 
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Equation  {3-112)  shows  that  Equation  (3-111)  is  satisfied,  therefore 
there  is  a solution.  This  solution  is 


K = [BpTBpJ-lBpT[A^  . . 


Again  the  solution  is  the  same  as  that  found  in  Examples  3.3  and  3.4. 
It  should  be  restated  that  there  is  not  always  a solution  to  the 
problem  via  Equation  (3-110)  and  to  insure  that  a solution  does  exist 
Equation  (3-111)  should  be  checked. 

This  process  can  be  extended  to  systems  where  the  input  matrix 
of  the  model  does  not  equal  that  of  the  plant. 


Proposition  3 
Given 


and 


Xp  = \Xp  + BpU 


X = A X + B u 
-m  m-m  m- 


where  Ap  and  A^  have  dimension  n x n.  Bp  and  B^  have  dimension  m x m, 
and  Bp  has  full  rank  then  there  is  a K matrix  and  a Q matrix  such  that 


[SI  - (Ap  t SpKlJ-'BpQ  = [si  - A„,r'B„ 


(3-114) 


if  and  only  if 


A,  + B [B  ^B  J-Ib  T[a  - a ] = a 
P PPP  pm  p m 


(3-115) 


74 


and 


Bp[Bp  BpJ  = Bn]  (3-116) 

If  Equations  (3-115)  and  (3-116)  are  satisfied  then, 

K = [BpV'®pX  ■ V 

and 

q - [Bp^Bp]-'8p\  (3-118) 

The  block  diagram  for  the  compensated  system  is  given  in  Figure  3-9. 


IV.  APPLICATION  OF  COMPENSATOR  DESIGN  ALGORITHMS  TO  THE 
AIRCRAFT  LONGITUDINAL  CONTROL  PROBLEM 


As  shown  in  Equation  (2-17)  the  aircraft  longitudinal  system  in 
state  space  form  has  an  A matrix  with  dimension  5x5  and  a B matrix 
with  dimension  5x4  where  the  matrices  are  given  in  Figures  2-2  and 
2-3  respectively.  Since  the  B matrix  has  two  rows  of  zeroes  it  has 
rank  less  than  or  equal  to  3 and  is  not  of  full  rank.  An  inspection  of 
the  B matrix  reveals  that  the  first,  second,  and  third  columns  are 
linearly  independent;  therefore  the  B matrix  has  rank  = 3.  Since  it  is 
necessary  for  the  input  matrix  to  have  full  rank,  as  shown  in  Propo- 
sition 2,  one  input  will  be  neglected  so  that  the  modified  B matrix 
will  have  full  rank.  After  the  compensator  is  designed,  the  neglected 
input  will  be  added  back  into  the  system.  Since  the  first  three 
columns  of  the  B matrix  are  linearly  independent,  U4  = aCq^  will  be 
neglected.  The  modified  system  now  is  of  the  form 

X = Ax  + Bu^ 

where 

A has  dimension  5x5 

B has  dimension  5x3 

The  X and  ^ vectors  are  now 
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r 


77 


AV 

Aa 


q 

aZ 


I 


I 

1 


(4-1) 


u 


(4-2) 


and  the  new  B matrix  will  be  the  first  three  columns  of  the  original  B 
matrix. 

If  the  system  is  controllable  with  controllability  index  p = 2, 
then  one  state  will  have  to  be  added  to  the  system  to  make  it  equi con- 
trollable (N=n+l=pm).  The  dimension  of  the  augmented  system  will  be 
N = 6. 


The  augmented  system  input  matrix  has  N-m  = 3 zero  rows,  two  from 
the  original  system  and  one  from  the  additional  state.  Theorem  I showed 
that  if  the  system  had  a controllability  index  p = 2 and  n-m  rows  of 
zeroes  in  the  B matrix,  then  the  similarity  transformation  matrix  T is 
dependent  only  on  the  B matrix  and  the  row  elements  of  the  homogeneous 
states.  Therefore  the  model  must  have  the  same  row  elements  as  the 
plant  in  the  third  and  fifth  rows  of  its  A matrix  and  the  additional 
state  must  be  added  in  the  same  manner  for  both  systems.  This  is  no 
problem  since  it  was  shown  in  Chapter  II  that  the  desired  response  could 
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be  obtained  by  adjusting  the  elements  of  the  first,  second,  and  fourth 
rows  of  the  A matrix. 

Thus  far  it  has  been  shown  that  the  modified  system 

1.  has  a B matrix  with  full  rank  = 3 and  n - m zero  rows. 

2.  can  be  altered  to  obtain  the  desired  response  without 
effecting  the  similarity  transformation  matrix. 

3.  will  have  to  be  augmented  with  one  additional  state  to  make 
it  equicontrol lable. 

Now  the  only  criterion  remaining  to  be  fulfilled  to  insure  a 
solution  to  the  problem  is  that  the  system  must  be  controllable  with 
controllability  index  p = 2.  Controllability  will  be  dependent  on  the 
values  of  the  A and  B matrix  elements  or  more  precisely  on  the  aircraft 
and  its  aerodynamic  coefficients. 

Examples  4.1  and  4.2  illustrate  the  algorithms  that  have  been 
presented  using  data  on  two  different  aircraft.  A computer  program 
was  written  to  compute  the  A and  B matrices  from  the  general  longitu- 
dinal aircraft  equations  and  the  aerodynamic  coefficients  of  the  air- 
craft in  question.  The  program  computes  the  system  transfer  function, 
the  roots  of  the  characteristic  equation,  and  the  short  period  and 
phugoid  damping  ratios  and  natural  frequencies.  The  program  was 
written  in  FORTRAN  IV  and  run  on  the  IBM-370  system  using  double  pre- 
cision arithmetic. 


f 
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Example  4.1 

Using  data  on  a business  jet  aircraft,  the  A and  B matrices  were 
determined.  They  are 


~7. 54(10)"^ 

8.47 

-32.2 

0 

0 

-2.08(10)'^ 

-.669 

0 

.996 

-4.96(10)"® 

0 

0 

0 

1 

0 

9.47(10)"^ 

-7.18 

0 

-1.35 

2.02(10)"® 

0 

675 

-675 

0 

0 

(4-3) 

”76.89 

0 

0 

-76.98" 

-5.36(10)”^ 

-.114 

0 

0 

0 

0 

0 

0 

(4-4) 

2.18(10)'^ 

.046 

11.64 

0 

0 

0 

0 

0 _ 

The  characteristic  equation  for  the  system  is 

C.E.  = + 2.02a^  + 8.07a^  + .0901A^  + .0935X  + 2. 57(10)"^  = 0 

(4-5) 


The  roots  of  the  characteristic  equation  are 


X.,  = -2.75(10)"^ 

X2,X3  = -2.77(10)"^  + i(.108) 
X^.Xg  = -1.01  + i(2.65) 


From  the  characteristic  roots  the  phugoid  and  short  period  damping 


80 


ratios  and  natural  frequencies  were  determined.  They  are 

= .026  , ojp  = .108  (4-6) 

P . 

= .355  , ojn  = 2.835  (4-7) 

sp  tigp 

and  L^/mV  were  computed  using 

L„/mV  = 1/2pV^SCl  /mV  (4-8) 

a 

and 

They  are 

IJmM  = .666 

Na  = 13.9 

Since  < 15,  the  performance  criteria  given  in  Figure  2-4  will  be 

used  to  evaluate  the  short  response  of  the  open  loop  plant  and  to 

determine  the  desired  u.  and  Figure  4-1  shows  the  location  of 

''sp 

the  response  characteristics  of  the  uncompensated  system  and  the  desired 
system  on  the  performance  criteria  plot. 


From  the  desired  location  on  the  performance  plot  given  in  Figure 

4-1,  (On  and  Ccn  the  model  can  be  determined  to  be 
"sp 
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Now  the  desirec  response  for  the  phugoid  mode  must  be  determined. 
From  Equation  (4-6)  is  noted  that  the  phugoid  period  is  approximately 
60  seconds  with  very  little  phygoid  damping.  From  the  performance 
criteria  it  was  shown  that  a minimum  damping  of  .2  is  desired.  In  the 
niodel  will  be  decreased  and  will  be  increased  to  approximately 
.4. 

To  develop  the  model  a22»  324’  ®42’  ®44  selected  in  the 

manner  given  in  Chapter  II.  Using  the  values  of  Ccn  and  a)„  given  in 

F "sp 

Equations  (4-10)  and  (4-11)  the  desired  short  period  characteristic 
equation  was  determined.  It  is 

C.E.sp  = + 2.1X  + 2.25  (4-12) 

First  a22  is  set  equal  to  a22  and  a24  is  set  equal  to  1.  Therefore 


ajjp  = .669 


Using  the  value  of  a22 

^m 

determined  by 


and  the  fact  that  a2^  ^ ^ ’ ®42  ®44 

m ^m 


344^  = + 2csp(On 


) 

sp 


.669  - 2.1  = - 1.43 


and 


342  = (322  )(a44  ) - ^ = .96  - 2.25  = - 1.29 

m m \ "sp 
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Adjusting  the  phugoid  roots  is  considerably  more  difficult  than  adjusting 
the  short  period  roots  as  shown  in  Chapter  II.  As  a general  rule  of 
thumb,  increasing  the  magnitude  of  the  n-ji  element  increases  the  phugoid 
damping  and  decreasing  the  magnitude  of  the  a^^*  ^■i3»  O’"  ^21  o^O'’’oot 
decreases  the  phugoid  natural  frequency.  The  phugoid  approximation 
given  in  Chapter  II  seems  to  break  down  as  the  phugoid  damping  becomes 
large. 


After  several  trys,  a satisfactory  model  was  determined.  The 
model  A matrix  is  given  in  Equation  (4-13). 


A 


m 


-.1 

8.47 

-20.0 

0 

0 

1 

0 

1 

-.669 

0 

1 

-4.96(10) 

0 

0 

0 

1 

0 

9.47(10)"^ 

-1.29 

0 

-1.43 

2.02(10) 

0 

675 

-675 

0 

0 

(4-13) 


The  model  characteristic  equation  is 

C.E.  = + 2.2X^  + 2.46X^  + .242x^  + .021x  + 5.6(10)'^  = 0 


and  its  pole  locations  are 
A*!  — — • 037 

A2.X3  = -.031  ± i(-.076) 
^4»^5  = “1*05  ± i(1.07) 
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and  the  damping  ratios  and  natural  frequencies  are 
Cp  = -38  , Wp  = .0825 

Csp  ~ '1  » “rigp  “ ^ 

was  made  the  same  as  Bp  to  insure  the  similarity  transformations  were 
the  same  betiveen  the  plant  and  model. 

Now  that  a satisfactory  model  has  been  constructed,  a compensator 
which  will  cause  the  aircraft  to  respond  in  the  same  manner  as  the 
model  will  be  designed.  First  the  compensator  will  be  designed  using 
Curran's  method.  The  compensator  was  designed  on  the  HP-2000  timeshare 
computer  using  Curran's  algorithm  as  outlined  in  Chapter  III. 

First  the  plant  was  augmented  and  placed  in  equicontrol lable 
canonical  form.  One  state  was  added  to  the  system  to  make  it  equicon- 
trollable  by  the  process  given  in  Proposition  1.  The  additional  state 
equation  is 


xg  = xi  - lOxg 

The  augmented  system  is  now  of  the  form 

A A A A 

X = Ax  + B^ 


where 


A 

A = 


I 0 
I 0 

A ' 0 

I 

I 0 
I 0 

0 0 0 0 [-10 


(4-14) 


(4-15) 


(4-16) 
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and 


B = 


B 

0 0 0 


(4-17) 


In  order  to  transform  the  system  to  equicontro] lable  canonical  form, 
transformation  matrix  Tp  was  generated.  It  is 


768.9 

0 

0 

76.9 

0 

0 

2.2(10)"^ 

.046 

11.64 

-5.4(10)“^ 

-.114 

0 

2.2{10)"^ 

.046 

11.64 

0 

0 

0 

0 

0 

0 

2.2(10)“^ 

.046 

11.64 

-3.62 

-76.9 

0 

0 

0 

0 

76.9 

0 

0 

0 

0 

0 

(4-18) 

The  plant  matrices  in  equi control lable  canonical  form  are 


0 

0 

0 

1 

0 

0 

0 

0 

0 

C 

1 

0 

0 

0 

0 

0 

0 

1 

-.076 

-.014 

-3.60 

-10.0 

-.013 

0 

1.42 

.270 

68.5 

.58 

-.66 

.419 

.055 

-.030 

-7.45 

.009 

.068 

-1.35 

(4-19) 

and 


AO-AO«7  Sl« 


UNCLASSIFIED 


AIR  FORCE  INST  OF  TECH  HRIOHT-FATTERSON  AFB  OHIO 
AIRCRAFT  COMPENSATOR  OESI0N  METHODS. (U) 

JUN  77  T L RIOOS 

AFIT-CI-77-I5 


F/®  20/A 


V 


NL 
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(4-20) 


Next  the  nx>de1  was  augmented  and  transformed  to  equicontrollable 
canonical  form.  The  state  was  added  to  the  model  in  the  same  manner 
as  the  plant  to  insure  the  transformation  matrices  would  be  equal. 

The  program  compares  the  similarity  transformations  for  equality  and 
exits  with  an  error  message  if  they  are  not  equal.  The  model  matrices 
in  equicontrollable  canonical  form  are  given  by  Equations  (4-21)  and 
(4-22). 


0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

-1.0 

o 

o 

1 

-1.75 

-10.1 

-.013 

0 

.735 

.269 

68.4 

.512 

GO 

1 

0 

.060 

-.006 

-1.56 

.006 

.010 

-1.43 

I 


(4-21) 


(4-22) 


\ 
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Using  the  method  given  In  Chapter  III  the  compensator  was  designed 
and  returned  to  the  original  coordinates.  The  compensator  output  was 
not  a function  of  the  additional  state  therefore  It  consisted  of  pure 
state  feedback.  The  feedback  matrix  Is 


-3 


-1.2(10) 

-8.93(10)' 


0 .159  0 

0 -7.47(10)'^  -.036 


-7.2(10)"^  0 


3.78(10)"®  .505  0 

(4-23) 

After  designing  the  compensator,  the  neglected  Input  Is  added 
back  Into  the  system  so  that  the  compensated  system  Is  of  the  form 
given  In  Equation  (4-24). 


X = [A  + BK']x  + Bu 


(4-24) 


where 


K'  = 


0 0 0 0 0 


and  A and  B are  the  original  system  matrices. 

The  feedback  matrix  was  also  determined  by  the  alternate  method 
given  In  Chapter  III  by  Proposition  3.  Using  the  matrices  given  In 
Equations  (4-3),  (4-4),  and  (4-13)  the  feedback  matrix  was  determined  by 


^ ■ [b/b  - A ] 


(4-25) 
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As  show  Proposition  3,  Equation  (4-25)  gives  a valid  solution  If 


and  only  r. 


- A„]  - A„ 

p P P , p ■■  m p-*  m 


(4-26) 


For  this  system 


B [B  ^B  ]"^B  ^ = 
P P*^  P 


1 0 0 0 0 
0 10  0 0 
0 0 0 0 0 
0 0 0 1 0 


0 0 0 0 0 


(4-27) 


From  Equation  (4-27)  It  Is  obvious  that  the  third  and  fifth  rows  of 
the  plant  and  model  must  be  equal  for  Equation  (4-26)  to  be  satisfied. 
This  Is  the  same  criterion  for  the  similarity  transformation  matrices 
to  be  equal. 

The  K matrix  using  the  alternate  method  Is 


-1.2(10)' 


K = -8.93(10)'^  0 -7.49(10)"^ 

3.78(10)"®  .505  0 


-7.2(10)"'’  0 


(4-28) 


Note  that  Equations  (4-28)  and  (4-23)  are  identical.  This  is  not 
surprising  since  the  solution  is  unique  regardless  of  the  method  used 
to  obtain  the  solution. 


89 


Using  CSMP  (Continuous  System  Modeling  Program),  the  uncompensated 

plant,  model,  and  compensated  plant  impulse  time  response  curves  were 

plotted  and  compared.  The  systems  were  perturbed  by  a delta  function 

on  inputs  aC.  and  ACm  at  t » .5  seconds.  Figures  4-2  through  4-11 
^c  c 

give  the  results  of  the  simultaneous  simulation.  Figures  4-2  through 
4-6  terminate  at  t = 15  seconds  to  compare  the  short  period  response 
of  the  system  states.  Figures  4-7  through  4-11  terminate  at  t = 250 
seconds  to  show  the  phugoid  response.  Note  that  the  model  response 
and  the  compensated  plant  response  are  the  same.  This  is  not  surprising 
since  the  purpose  of  the  feedback  matrix  in  the  compensated  system  is 
to  equate  the  model  and  compensated  system  transfer  functions. 

Example  4.2  illustrates  that  once  a compensator  is  designed  which 
will  cause  the  aircraft  response  to  be  well  within  the  satisfactory 
regions  of  the  performance  criteria  at  a median  velocity  and  altitude, 
then  it  will  remain  within  that  satisfactory  region  over  a range  of 
flight  conditions.  Wind  tunnel  data  of  the  various  aerodynamic 
coefficients  for  a jet  fighter  type  aircraft  were  used  to  construct  the 
system  matrices  over  a range  of  velocities  and  altitudes.  From  the 
aerodynamic  data,  interpolating  polynomials  that  were  a function  of 
velocity  and  altitude  were  written  for  each  of  the  aerodynamic  coeffi- 
cients. The  interpolating  polynomials  were  accurate  for  velocities 
ranging  from  330  fps  to  880  fps  and  altitudes  ranging  from  sea  level 
to  30,000  ft.  The  system  was  simulated  only  in  these  ranges  of  flight 
conditions.  A modified  version  of  the  program  used  in  Example  4.1 
was  used  for  this  program.  Using  the  interpolating  polynomials,  the 
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Figure  4-2.  Change  in  Velocity  vs.  Time  For  System  in  Example  4.1. 


Figure  4-3.  Change  in  Angle  of  Attack  vs.  Time  For  System  in  Example  4.1 


Figure  4-4.  Change  in  Pitch  Angle  vs.  Time  For  System  in  Example  4.1. 


Figure  4-5.  Pitch  Rate  vs.  Time  For  System  in  Example  4.1 
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Figure  4-6.  Change  in  Altitude  vs.  Time  For  System  in  Example  4.1. 
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Figure  4-7.  Change  In  Velocity  vs.  Time  For  System  in  Example  4.1. 


Figure  4-8.  Change  in  Angle  of  Attack  vs.  Time  For  System  in  Example  4.1. 


Figure  4-9.  Change  in  Pitch  Angle  vs.  Time  For  System  in  Example  4.1. 


Figure  4-10.  Pitch  Rate  vs.  Time  For  System  in  Example  4.1. 
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program  calculated  the  aerodynamic  coefficients  and  then  the  system 
matrices  for  a range  of  velocities  and  altitudes.  The  program  tabu- 
lated N^,  (u^p,  and  the  phugoid  roots  for  velocities  ranging  from 

330  fps  to  880  fps  in  110  fps  increments  and  altitudes  from  0 feet  to 
30,000  ft.  in  5000  ft.  increments  at  each  velocity. 


Example  4.2 

The  velocity  and  altitude  at  which  this  aircraft  compensator  was 
designed  was  V = 660  fps  and  Z = 0 ft.  The  system  matrices  at  these 
flight  conditions  are 


1.45(10)‘ 


5.75(10)' 


16.45 

-1.19 


-8.10 


-32.15 


■1.26(10)' 


5(10)' 


(4-29) 


257.94 


-259.03 


-.036 


34.29 


The  characteristic  equation  for  the  system  is 


(4-30) 


C.E.  = \^  + 2.43X^  + 9.50X^  + .156X^  .047X  1.14(10)"^  = 0 


(4-31) 
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and  the  roots  are  located  at 
= -2.44(10)’^ 

= -6.38(10"^)  + i(.07) 
x^.Ag  = -1.21  + i(2.83) 

The  short  period  and  phugoid  damping  ratios  and  natural  frequencies 
are 


and 


Cgp  = -39  , 0)^  =3.08 

= .09  w = .07 

P "n 


N = 24.27 
a 

Figure  4-12  shows  the  intersection  of  and  i'or  the 

uncompensated  plant  and  the  proposed  compensated  plant. 

From  Figure  4-12  it  can  be  seen  that  the  model  short  period  ^ 
and  0)^  will  have  to  be 


C 


sp 


.7 


(4-32) 


0)^  = 2.427  (4-33) 


The  phugoid  damping  for  the  model  must  also  be  increased  for  a satis- 
factory response.  Again  a phugoid  damping  of  approximately  .4  will  be 
satisfactory. 
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The  model  was  designed  in  the  same  manner  as  the  model  in  Example 
4.1.  The  model  A matrix  is 


-.0625 

3.30 

-32.15 

0 

0 

-1.83(10)"^ 

-1.19 

0 

1 

-1.26(10) 

0 

0 

0 

1 

0 

5.75(10)"^ 

-3.27 

0 

-2.21 

5(10) 

0 

660 

-660 

0 

0 

(4-34) 


The  model  characteristic  equation  is 

C.E.  = + 3.46X^  + e.lOX^  + .373X^  + .025X  + 1. 95(10)"^  = 0 

and  the  model  characteristic  roots  are 


X^  = -8.91(10)"^ 

^2*^3  ~ + i(.055) 

X^,Xg  = -1.7  ^ i(1.73) 

The  model  short  period  and  phugoid  damping  ratios  and  natural  frequen 
cies  are 


^sp  ■ 

.7  . 

u)„  = 2.427 

"sp 

.43 

u = .061 

sp 
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The  feedback  matrix  was  determined  by  both  Curran's  method  and  the 
alternate  method  with  the  same  result.  The  solution  is 


”-1.80(10)"^ 

-.051 

0 

0 

6" 

1.14(10)’^ 

4.66(10)"^ 

0 

-.032 

0 

(4-35) 

-4.44(10)'^ 

.141 

0 

-.029 

0_ 

The  uncompensated,  model,  and  compensated  system  were  simulated 
using  CSMP  as  in  Example  4.1.  Figures  4-13  through  4-22  give  the 
results  of  the  simulation.  Again  note  that  the  compensated  system  and 
the  model  time  response  are  identical,  indicating  their  transfer 
functions  are  equal. 

From  the  tabulated  values  of  N^,  , and  the  compensated 

aircraft  peformance  was  plotted  on  graphs  corresponding  to  Shomber 
and  Gertsen's  performance  criteria.  Figures  4-23  through  4-30  show 
the  results.  Each  plot  is  at  a fixed  velocity  while  the  altitude  is 
varied.  The  plots  indicate  that  the  aircraft  performs  satisfactorily 
for  all  flight  conditions  plotted. 
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Figure  4-20.  Change  in  Pitch  Angle  vs.  Time  For  System  in  Example  4.2 


Figure  4-22.  Change  in  Altitude  vs.  Time  For  System  in  Example  4.2. 


Figure  4-23.  Handling  Qualities  of  Compensated  System  In 

Example  4.2  for  Velocity  = 330  fps  and  Altitude 
Ranging  From  Sea  Level  to  30,000  ft. 
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Figure  4-24 


Handling  Qualities  of  Compensated  System  In 
Example  4.2  for  Velocity  = 880  fps  and  Altitude 
Ranging  From  Sea  Level  to  30,000  ft. 
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Figure  4-25.  Handling  Qualities  of  Compensated  System  In 

Example  4.2  for  Velocity  =440  fps  and  Altitude 
Ranging  From  Sea  Level  To  30,000  ft. 


Unaccertabl e 


•Acceptable 


Satisfactory 


Handling  Qualities  of  Compensated  System  In 
Example  4.2  For  Velocity  = 770  fps  and  Altitude 
Ranging  From  Sea  Level  To  30,000  ft. 
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Figure  4-27.  Handling  Qualities  of  Compensated  System  In 

Example  4.2  for  Velocity  = 550  fps  and  Altitude 
Ranging  From  10,000  ft.  to  30,000  ft. 


Unacceptabl e 


Acceptable 


Figure  4-28.  Handling  Qualities  of  Compensated  System  in 

Example  4.2  For  Velocity  = 550  fps  and  Altitude 
Ranging  From  Sea  Level  to  5,000  ft. 
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Figure  4-29.  Handling  Qualities  of  Compensated  System  in 
Example  4.2  for  Velocity  = 660  fps  and  Alti 
Ranging  From  20,000  ft.  to  30,000  ft. 


Handling  Qualities  of  Compensated  System  in 
Example  4.2  for  Velocity  = 660  fps  and  Altitude 
Ranging  from  Sea  Level  to  15,000  t. 
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V.  SUMMARY  AND  CONCLUSIONS 


The  longitudinal  aircraft  equations  were  investigated  to  determine 
which  elements  of  the  A matrix  dominate  the  response  of  the  aircraft. 

A method  for  contriving  a model  that  has  a desired  response  based  on 
performance  criteria  was  developed.  Two  methods  for  designing  a 
compensator  which  forces  the  aircraft  to  respond  in  the  same  manner  as 
the  desired  model  were  presented. 

The  method  for  designing  the  model  is  straightforward.  The 
system's  short  period  response  is  easily  and  exactly  adjusted  by  altering 
three  elements  of  the  A matrix.  The  phugoid  response  was  harder  t 
adjust  precisely  because  of  the  indefiniteness  of  the  dominant  phugoid 
elements  in  the  A matrix.  This  resulted  in  a trial  and  error  method 
for  obtaining  the  exact  desired  phugoid  response.  However,  it  should 
be  noted  that  this  difficulty  in  adjusting  the  phugoid  response  is  not 
a major  problem  since  the  phugoid  criterion  is  not  precise. 

A special  case  based  on  equality  of  the  plant  and  model  similarity 
transformations  in  Curran's  algorithm  was  presented.  The  criteria  for 
equality  of  the  similarity  transformations  were  derived  and  it  was 
shown  that  the  aircraft  equations  in  state  space  form  met  the  criteria 
if  the  systems  are  controllable.  The  advantage  of  Curran's  algorithm 
over  other  methods  using  optimal  control  theory  is  the  directness  of 
solution.  The  algorithm  was  easily  implemented  in  a computer  program 
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and  a solution  to  the  program  involved  a minimum  of  human/computer 
interaction. 

The  second  method  for  designing  the  compensator  was  more  direct 
than  Curran's  method  since  it  did  not  employ  any  similarity  transforma- 
tions nor  was  it  necessary  to  augment  the  system  as  in  Curran's  method. 
The  alternate  method  imposed  the  same  restrictions  on  the  system 
matrices  as  Curran's  method,  e.g.  the  B matrix  must  have  full  rank  and 
the  model  must  be  structurally  similar  to  the  plant.  It  was  shown  that 
a necessary  condition  for  the  similarity  transformations  to  be  equal 
in  Curran's  method  was  that  the  model  and  plant  input  matrices  had  to 
be  equal.  The  alternate  method  did  not  impose  this  restriction. 

There  were  restrictions  on  the  alternate  method  involving  the  structure 
of  the  model  matrices  and  t plant  input  matrix.  They  are  presented 
in  Proposition  3 in  Chapter  III. 

The  restriction  that  the  input  matrix  have  full  rank  was  shown 
to  be  no  limitation  for  the  compensator  design  for  any  controllable 
system.  In  Proposition  2,  it  was  shown  that  the  input  matrix  can  be 
made  full  rank  by  not  utilizing  all  of  the  control  inputs. 

In  Chapter  IV,  compensators  were  designed  for  two  aircraft  using 
both  Curran's  method  and  the  alternate  method.  The  results  were  more 
than  satisfactory.  Simulation  results  verified  that  the  compensated 
system  response  was  equal  to  the  model  response.  It  was  also  shown 
in  Example  4.2  that,  once  the  aircraft  is  compensated  so  that  its 
response  at  some  arbitrary  median  flight  condition  is  well  within  the 
satisfactory  region  of  the  performance  criteria  plots,  then  the  system 
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response  will  be  satisfactory  over  a range  of  flight  conditions. 

Example  4.2  verified  the  derivation  given  in  Equations  (2-34)  through 
(2-43)  for  determining  how  the  compensated  system  and  varied 
with  velocity  and  altitude.  The  fact  that  the  compensated  system 
response  is  satisfactory  for  a range  of  flight  conditions  suggests  that, 
for  aircraft  that  have  a large  flight  envelope  such  as  supersonic  or 
high  altitude  aircraft,  several  compensators  could  be  designed,  each 
for  overlapping  flight  regimes.  This  would  involve  some  sort  of 
sensing  and  switching  control  circuitry  to  select  the  proper  compensator 
for  the  corresponding  flight  regime. 
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